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Abstract We use Renewal Theory for the estimation and interpretation of the 
flare rate from the Geostationary Operational Environmental Satellite (GOES) 
soft X-ray flare catalogue. It is found, that in addition to the flare rate variability 
with the solar cycles, a much faster variation occurs. The fast variation on time 
scales of days and hours down to minute scale appears to be comparable with 
time intervals between two successive flares (waiting times). The detected fast 
non-stationarity of the flaring rate is discussed in the framework of the previously 
published stochastic models of the waiting time dynamics. 

Keywords: Solar Flares, Waiting Times, Flaring Rate, Reliability Theory, Haz- 
ard Function 

1. INTRODUCTION 

The phenomena related to energy transformation and release on the Sun, possi- 
bly are of highest importance for modern solar astrophysics. Solar flaring, as a 
steady process of energy release, plays a central role and has been drawing the 
attention of the scientific community for decades. 

In this work, we focus on the statistical properties of so-called solar flares wait- 
ing time(s) (wt), i.e. the interval between two flares close in time. The data are 
provided by the Geostationary Observational Environmental Satellites (GOES) 
in the soft X-ray band. This catalogue is chosen as the longest available record 
of uninterrupted observations. The Ha flare records from the NGDC-NOAA 
catalogu^] are shorter, and used to emphasise the invariance of the reported 
results. 
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Flare waiting times statistics has been extensively debated in the literature. 
Boffetta et al. (19991 fitted the waiting time probability density function (pdf^ 



with a power-law in the range between 6 and 67 hours. This estimate was im- 
plemented using a single record over 20 years. Wheatland ( 2000[ ) considered the 
variation of the mean flaring rate with the solar cycle, and proposed the model 
of the time-dependent Poisson process (see also Moon et al. 2001 ). In turn, this 



model gives a power-law-like pdf only in the limit (Wheatland and Litvinenko, 



2002). However, Lepreti, Carbone, Veltri (20011 argued for a local departure of 
the wt series from the Poisson process, i.e. "memory" had been detected in the 
data. 

We apply an alternative statistical description to the methods that had been 
used in the articles cited above. This method is equivalently unambiguous in 
describing the stochastic processes. The flare rate is estimated explicitly from 
the waiting time pdf. Such an approach has the significant advantage of linking 
the waiting time of the flare and its rate, which is the instantaneous probability 
of the flare per unit of time. 

The paper is organised as follows. In Section [2] the rather simple mathemat- 
ics used in the paper is summarised, considered for the paper's consistency to 
avoid referring the reader to specific literature. The main results are presented 
Section [3] and interpreted in Section [4j The conclusion and final remarks are 
presented in Section [5] 



2. THEORETICAL BACKGROUND 



We adopted the terminology of Renewal Theory (e.g. Cox 1962) to set a math- 
ematical framework and for intuitive and easy further interpretation. 

We focus on two objects of study: random events and time interval between 
near-in-time events, i.e., waiting times. By analogy with Renewal Theory, the 
failure of some abstract device is attributed as the elementary random event, 
which we assign to the flare. Then the flare waiting time is associated with a 
random non- negative variable X, called failure tim^\ being the interval between 
adjacent]^] failu res, that is, between two flares. 

The random variable X is characterised by a probability density function f(x) 

,. . Prob(x < X < x + Ax) 

f(x) = lim , (1) 

J y 1 Ai->o Ax v ' 

with f(x)dx = 1. 

The probability that a flare has occurred (a device has failed) by time x is 
given by the cumulative distribution function F(x): 



2 We use the mathematically strict definition of the density as the differential quantity with 
respect to the distribution, which is often used in a misleading way instead. 

3 This time interval is also called the age of a device, i.e., its lifetime without any failure. 

4 An abstract picture is considered here: a broken device replacement (repair) is thought to 
be a very instant, in such a way that the time required for the repair is zero. A continuous 
operation is meant, and provided, for instance, by multiple hardware availability. 
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F(x) = Prob(X < x) = I f(u)du. (2) 



o 



The probability that a flare has not occurred (a device has not failed) up to 
time x is given by the survivor function J°(x) : 



F(x) = Prob(X >x)= (3) 
l-F(x) = / f{u)du, (4) 



The probability of immediate failure of a device (occurrence of a flare event) 
known to be of age x (no flare during x) is given by the age-specific failure rate 
h(x^ Consider a device known not to have failed at time x and let h(x) be the 
limit of the ratio to Ax of the probability of failure in time interval (x, x + Ax] 

Prob(x < X < x + Ax\x < X) 

nix) = lim t . (5) 

v ; a^o Ax 

The latter permits further transformation according to the definition of the 
conditional probability for two events a, b, 

„ w ,,n Probia and b) 

Pr ° b ^= Prob(b) - (6) 

The event 

"x < X < x + Ax and x < X" 
is included in essentially the same event 

"x < X < x + Ax" ; 

then Equation ^ reads 
h(x) = 

Prob(x < X < x + Ax) 1 



= lim 

Ano Ax Prob(x < X) 

= Six) 
F{x) 

and the next form is used for the computations in the paper: 



5 This function is also known as the hazard function, hazard rate or failure rate. 
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h ( x ) = M = IM ( 7 ) 

The stochastic Poisson process has the exponential waiting time pdf 

/(*) = Xe~ Xx , (8) 
where A =< x > is the constant failure rate 

M*) = ^ Xudu = A- 0) 

Jx 

Conversely, if the stochastic process has a constant h, it is a Poisson pro- 
cess. The constancy of h p (x) reveals the "no-memory" property of the Poisson 
stochastic process. 



3. DATA ANALYSIS 

The GOES flare catalogues are analysed with particular interest in the interval 
between consecutive flares. The GOES observations provide the longest record 
and complementary Ha data provided allow one to consider both X-ray and 
Ha flares. 

The catalogued event is recorded by specifying starting time, peak time and 
ending time. Thus, the waiting time is not unambiguously defined. We consider 
two definitions: the Peak-to-Peak (PtP) waiting time — the interval between 
times of maximum flux rate of two near-in-time flares, and the End-to-Start 
(EtS) waiting time — the interval between the end time of the predecessor and 
the starting time of the current event. The latter definition is intuitive, but might 
be of questionable applicability from the point of view of statistical methods that 
will be used (see Section [4]). However, it is used in the article to demonstrate the 
influence of the definition of the wt on its statistical properties. 

Initially, the catalogue should be preprocessed to avoid errors due to data 
gaps. The neighbourhood events to every catalogue item marked with the data 
gap flags (notably D, E) were ignored; thus, after filtering we have the set of 
52181 X-ray events (including GOES B-class events) dating from September 
1975 until December 2008. The Ha events record is shorter, viz. 11124 events 
and it is used only for numerical comparison in the following. 

The calculated PtP waiting times from the GOES catalogue are shown in Fig- 
ure[l] The correlations between monthly averaged sunspots number and monthly 
averaged wt in both both Soft X-rays and Ha are evident and intuitive. During 
the solar minimum years the flare rate is lower, leading to longer wt prevalence; 
conversely, during the solar maximum years the rate is higher and thus waiting 
times are shorter on average. 

This substantial variation leads to separate consideration of the solar cycle 
phases according to corresponding traces in the calculated wt series. To separate 
phases, we consider the time derivative of the monthly average sunspot number 
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N = dN(t)/dt. The qualitative agreement (in Figure[l]is shown by grey lines) in 
the fluctuations of wt and N(t) is attained empirically, by considering intervals 
with limited fluctuations in N around the corresponding zeros, to be less than 
25% for the minima phases and less than 50% for the maxima. 

Figures 2][4 show a pdf of waiting times for joint solar cycle phases, solar min- 



imum and maximum, respectively, by means of the two wt definitions mentioned 
above. The solid line in each panel stands for the reference exponential distribu- 
tion \e~ Xx , with the average A" 1 estimated from the series whose pdf is shown 
by the scatter plot. To some extent, the similarity with the pure exponential is 
attained in a single case only (the bottom panel in Figure |4| . 

The density shapes are similar within a given wt definition. Remarkably, the 
variation in wt definition affects the most probable events: the bell-like pdf versus 
plateau- like over almost 1.5 decades. 

The plots in Figure [2] should be interpreted as power-laws according to |Bof-| 



fetta el al. (1999). However, the power-law fit has to be considered with special 
care, as the findings illustrated in the following are challenging the straightfor- 
ward applicability of the power-law fit. 

3.1. Flaring rate estimation 



Wheatland (20001 had ignored flares whose peak X-ray flux is less than 10 -6 
Wm" 2 , due to substantial variation of the soft X-ray background with the solar 
cycles. Thus GOES .B-class events (f?-flares for short) had not been considered. 



Following Wheatland (2000), additional sets of data without _B-flares were 
generated. In the following, the effect of this removal is considered systematically. 

We estimate the flaring rate h explicitly from the pdf according to Equa- 
tion ([7]). The result for joint solar cycle phases is reported in Figure [5j in Figure [6] 
joint phases with £?-flares excluded are shown for comparison, and the phase- 
wise estimated rates are shown in Figure [7j The effect of the wt definition on 
the curvature of the estimated function h is shown for the entire dataset. 

In Figure [5] the solid lines represent the smoothing adjacent running average 
with 10-point window, to emphasise the character of the functions. 

Figure [6] is complementary to the previous one, and demonstrates the alter- 
ation of the smoothed rates, when the S-flares are removed. The smoothed rates 
are shown by the scatter plot. 

In Figure [7J the estimated rates are shown for the different solar cycle phases 
separately. The estimates are smoothed by the 15-point adjacent running aver- 
age. The phases are coded by the colour of the scatter plot; the wt definition by 
the symbol shapes; the cases without i?-flares are shown by lines. 

The estimated flare rate allows one to reconsider the applicability of the 
power-law fit for the waiting time pdf by trivial algebra: the simplest functional 
form for the power-law fit is given by 

g[x)=Ax~ a , (10) 

where A and a are the estimated parameters, and due to Equation Q the 
corresponding failure rate to the power-law pdf g(x) reads 
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Table 1. Numerical characteristics of the datasets generated according to the 
PtP flare waiting time definition. 

Catalogue Phase < x > a [h] < x > b [h] -y a 7 b 



joint 


4.48 


5.89 


-0.86 ±0.006 


-0.91 ± 0.005 c 


X-ray max 


3.14 


3.27 


-0.65 ±0.008 


-0.73 ± 0.007 c 


min 


9.58 


28.99 


-0.61 ± 0.01 


-0.67 ±0.01 


joint 


6.00 


9.87 


-0.54 ± 0.006 


-0.53 ± 0.008 


Ha max 


4.58 


7.01 


-0.46 ±0.01 


-0.43 ±0.01 


min 


11.15 


50.14 


_d 




a events including GOES B-class flares. 

k events excluding GOES S-class flares. 

c the power-law region is remarkably present. 

rf insufficient record length to produce a reliable result. 






Table 2. Numerical characteristics of the datasets generated according to the 
EtS flare waiting time definition. Note: the power-law regions in these data appear to 
be more pronounced and longer. 


Catalogue Phase 


< x > a [h] 


< x > b [h] 




7 b 


joint 


4.15 


5.54 


-0.82 ±0.005 


-0.99 ± 0.006 c - d 


X-ray max 


2.79 


2.92 


-0.65 ±0.01 


-0.78±0.009 c 


min 


9.33 


28.72 


-0.66 ±0.01 


-0.56 ±0.007 


joint 


5.68 


9.54 


-0.50 ± 0.005 


-0.52 ±0.006 


Ha max 


4.20 


6.68 


-0.46 ±0.01 


-0.45 ±0.01 


min 


10.94 


49.75 


_e 


_e 



events including GOES S-class flares, 
events excluding GOES B-class flares, 
the power-law region is remarkably present. 

variations of the fit parameters may lead to the —1 power-law index, 
insufficient record length to produce a reliable result. 



Ax~ a x~ a 

with |a| > 1. 

To verify this relation for the estimated rates, we fit the smoothecj^] rate h(x) 



by a function of the form of Equation ( 10 1; the estimated exponent 7 is compared 
with —1. 

Tables [T] and [2] report the average wt and the exponent 7 for all datasets that 
had been analysed in this work. 



6 The smoothing is the same as for the graphs in 



Figure [7] 
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4. DISCUSSION 



In this section, the opposing results and interpretations existing in the literature 
are compared and discussed in the context of the findings reported. We shall 
highlight the discrepancies and agreement with those previously published by 
the cited authors. 

4.1. Uncertainties 

Evidently, our results are strongly affected by the completeness and accuracy 
of the GOES catalogue. In fact, the event detection is done in an automated 
way and any change in the detection algorithm would alter the statistics. Fur- 



thermore, the flare event obscuration (for details see |Wheatland 2001 ) brings 



about quite hard arguments against the reliability of any study: the missing 
of substantial amount of events (up to 75 ± 23% for events above GOES Cl- 
class) might be so crucial that the pdf would dramatically change its behaviour. 
However, our results are comparable with those based on the same datasets. 

In the case of the solar minima, the longer record is required most of all. 
The relatively high probability of longer wt (Figure [3]) conceals the probability 
variation along the domain, since the pdf is shrunk along the ordinate and has 
too wide spreading of the tail. The longer series should clarify the pdf shape. 

Equation ^ diverges as the cumulative distribution approaches values at the 
very tail 

lim F(x) = 1 , (12) 

x— J-oo 

and, consequently, the wide tails of the pdf are amplified by means of the non- 
linear transformation given by Equation Q. Thus, the resulting tails of the 
estimated rates appear to be even more scattered (Figure [5]) . This motivated 
us to apply a smoothing to the functions obtained, which, in turn, eliminates 
details and provides rather qualitative result. 

4.2. Waiting time definition 

A random evenlQ should be defined as a point-like instance in time i.e., the event 
has no length. From this point of view, the PtP definition is more mathematically 
adequate, rather than the EtS waiting time defined by the subtraction of the 
flare duration from the base time line. Practically, the flux maximum point is 
easier to detect precisely, with respect to the starting/ending times, since their 
flux values are closer to the background noise level. 



Wheatland, Sturrock, McTiernan (1998) and Boffetta et al. (1999) used the 



wt defined as the difference between the times of peak flux rates of two near- 



in-time flares. Later in Wheatland (20001 the wt was defined as the difference 



between the start times. However, the latter definition is qualitatively coincident 
with what we called PtP . We would ignore possible numerical discrepancies. 



7 Here by "event" we mean the event as denned in Probability Theory. 
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The major differences in the estimated rates caused by the wt definition had 
appeared in the range of short waiting times (Figures 5 — 7). 

In spite of the cut of wide-spreading tails, according to the authors the con- 
vergence of rates at longer waiting times can be identified, regardless of the 
wt definitions that had been used. This is seen where the tails of the flare rates 
seem to coincide. Such a convergent behaviour is independent of the solar cycle 
phase, as well as of the i?-flares exclusion (Figure [7]). However, the B-flares 
exclusion leads to more ambiguous plots, since fewer events are considered in 
the dataset. The joint consideration of the phases demonstrates this behaviour 
too (Figure [5]). 

If we assume that energetically large flares are separated by longer waiting 
times on the average, we can conclude that the relatively long duration of en- 
ergetic flares appears to be statistically indistinguishable from the long waiting 
times defined by the point-like events. Thus, the EtS definition tends to mimic 
the PtP in the right-hand half of the x domain, where the waiting times are 
longer. 

4.3. Waiting time pdf versus flaring rate 

In this work, the pdf are presented mostly for illustrative purposes, i.e. to 
demonstrate a degree of divergence between the exponentials, the variation due 
to the solar cycle phase and the effect of the wt definition. We do not consider 
the pdf fit, but study the underlying models by means of the estimated flaring 
rate. 

The failure rate formalism described in Section [2] leads to a broader view 
on the solar flaring and physics underlying the waiting time statistics. The 
function h(x) permits a qualitative description of the dynamics of the stochastic 
process being characterised by the pdf f(x). In fact, the rate h{x) and the pdf 
f(x) are related to different stochastic processes: the rate h(x) is the almost 
instantaneous probability of an event]^] with no event during the waiting time 
x, by the definition. In turn, the pdf f{x) is the probability change rate of the 
waiting time x. The function h(x), by definition, explicitly combines the almost 
instantaneous probability of the flare with the waiting time that had elapsed 
before its occurrence, i.e. it is the conditional probability. 

In other words, the rate h{x) gives the probability of the event delimiting the 
length of x, namely it is related to the process inducing the one being presented 
by the series of the waiting times in the catalogue. 

In the present work, we use the estimated rates h to analyse the flaring dynam- 
ics which can be derived from the GOES catalogue, and then we reconsider the 
cited stochastic models of the waiting time pdf fit. Particularly, in the piece- wise 
constant Poisson model the flare rate is a free parameter. 



per unit time (minute). 
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4.4. Non-constancy over the waiting time domain 



The most important feature of the estimated h is the variation over the x domain. 
The function h is the explicit non-linear function of x: 



with bell-like shape in semi-logarithmical plot. This has far-reaching conse- 
quences for the flaring dynamics. 

First, the flaring rate is by no means stationary: it non- linearly depends on 
waiting time of short scales of days, hours and minutes (the fast non-stationarity 
for short). 

Second, the flare rate has a certain behaviour: less energetic flares (£?-flares 
at large, but not only these) exhibit an increasing rate. Next, the rate reaches 
maximum values, and then it rather slowly declines (similarly to a power-law, 
which in some cases manifests itself very notably). Thus, we can point out char- 
acteristic waiting times, which indicate the time intervals of the most probable 
flare occurrence. We would not define this time as numerically accurate due to 
significant uncertainties. However, one can empirically define a range just by 
examining Figure [6] and Figure [7j for joint phases, say 35-45 rain. (EtS ); 50- 
60 min. (PtP ); for the solar minimum, say 20-55 min. (EtS ); 50-75 min. 
(PtP ); and for the solar maximum 30-45 min. (EtS ); 60-100 min. (PtP ). 

The character of the fast non-stationarity is invariable within the solar cycle 
phases (Figure [7]), and remains recognisable for smaller data records, when the 
B-flares are excluded (with a word of caution on the reliability of the result 
based on the smaller number of samples). 

4.5. Rate change with the solar cycle phase 

The slow non-stationarity of the flare rate is given by the solar cycles, say "large 
scale". The phases have different probability of flaring per unit time, as it is 
natural to assume just on the basis of the comparative overview in Figure [T] 

The most frequent and thus very short flares (corresponding to the smaller 
a;) are common features for both phases. They can hide statistical differences 
between solar cycle phases: this fact is supported by the match in the increas- 
ing regions of the rates corresponding to the disjoint phases. In Figure [7J the 
PtP rates are equivalent in the range from 4 to 20 min., regardless of the solar 
cycle phase considered. The same statement but somewhat weaker holds for the 
EtS rates. 

Qualitatively, the rate shape is invariable during the solar cycle, however, the 
variation in the mean wt may reach one order of magnitude. This is the case 
when the B-flares are excluded (see Tables [I] and [2]) . 

Summarising, the flare rate takes the form 



with explicit dependence on waiting time and generic time, i.e. the time frame 
when the corresponding observations were made. 



h = h(x) 



(13) 



h = h(x, t) 



(14) 
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4.6. GOES B-class events 



We considered separately the case of the excluded £?-flares. For joint phases 
almost 32% are B-flares and 3.4% for the solar maxima. 

Qualitatively, the rates retain their shape over the wt domain. Nevertheless, 
the rates without B- flares are systematically lower in the joint phases (Figure [6|) 
and with almost negligible variations during the solar maxima (in FigureT7^ 
dashed and dash-dotted lines). 

Substantial variation occurs during the cycle minima phases, when 80% of 
the events are B-flares. In this case, the notable rate matching in the shorter 
increasing range is very small for the PtP rate, and does not takes place for the 
EtS rates (in Figure [7] solid grey lines) . 

In fact, the elimination of the frequent events from the dataset is a very 
significant operation. Particularly, the power-law exponents are sensitive to the 
relative strength of the frequent events. This is caused by the intrinsic divergence 
of the power-law statistics at the origin, at the most probable {i.e. the most 
frequent) events. 

In some cases, excluding f?-flares decreases the reliability of the results be- 
cause of the small size of the resulting dataset. Nevertheless, it underlines the 
average waiting time order of magnitude change, which is the reason for the 
separate considerations of the solar cycle phases (see Tables [T] and [2] for numerical 
estimates). 

4.7. Power-law fit 



Considering Equation (11), we compare the flare rate power-law fit with the 
value of —1. 

One should notice that the tails of h{x) are quite uncertain in considering the 
indices listed in Tables [T] and [2] to be steady and exact: the change of the region 
boundaries being chosen for the fit can substantially modify the numerical value 
of the exponent, and, in a certain sense, it may be considered to be subjective. 

But we rely on the fact that the modulus of the estimated indices is systemat- 



ically less than 1. Thus, arguments for fitting the GOES data by Equation (10) 
can, very likely, be rejected. 

On the other hand, jointly considered solar cycle phases in the X-ray band 
with flares above B-class have revealed a value comparable with — 1. By noting, 
that this takes place when average waiting times differ by one order of magnitude, 
we would point the accordance with the power-law limit pdf of the non-stationary 
exponential random variable, reported by Wheatland and Litvinenko (2002). 



This argument is missing in Boffetta et al. (1999), when joint phases have been 



fitted confidently by the power-law. 
4.8. Excess of short events 

The excess of short waiting times was pointed out by |Pearce, Rowe, Yeung| 



( 1993) and Wheatland, Sturrock, McTiernan ( 1998 1. It appears to be an intrinsic 



property of the flaring process, and it is detectable regardless of the instrument 
and(or) the band that is used. 
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The rising character of the rate h at small x appears to be an indicator of the 
importance of the energetically small events. Short waiting times mostly separate 
less energetic events. Empirically, this point is supported by the influence of 
the 5-flares' presence on the flare rate. In turn, the increasing probability at 
short wt suggests that the true flare occurrence per short time unit is very 
large, but it is hidden due to sensitivity of the instruments and the obscuration. 
Certainly, somewhere on these time scales the very sympathetic flaring ( |Moon e\ 
al., 2002 Biesecker, Thompson, 2000| would take place, and one can realise that 



a dependence of the form h = h(x) poses explicitly the "memory" in the re cord^] 
Thus, it is hard to ignore the arguments by |Lepreti, Carbone, Veltri ( 2001 ) 
concerning the memory effect in the waiting time pdf. That is also in accordance 
with the divergence (viz. short flare overabundance) of the flare duration pdf at 
the origin, which is commonly accepted to be a power-law. 

The GOES soft X-ray flare catalogue gives an overall flaring picture, consid- 
ering the Sun as a whole. The beauty and power of the non-stationary Poisson 
model is in its potential ability to represent the global solar flare dynamics as 
consisting of somewhat trivial ("memoryless") mathematical objects — expo- 
nentials. Presumably, the different active regions (or even smaller regions) may 
contribute to the global flaring with notably different rates and thus appear in 
the limit as a power-law-like distribution. However, the individual active regions 
have very poor statistics (about 100 events at most) and a series of assumptions 
should be made a priori ( |Wheatland| |2001[ ). In addition, |Wheatland| fl2001[ ) 
reported a piece-wise constant Poisson fit for the active region, which gives a 
motivation to consider even smaller flaring areas to be an "elementary piece" in 
the mechanism just speculated on above. 



5. CONCLUSIONS 



The fast variation on time scales from minutes to hours of the flare rate is 
a highlight among other findings in this work. The large-scale variation with 
the solar cycle exhibits a complexity that can be coped with by a phase-wise 
splitting of the data. However, the fast non-stationarity during the waiting time 
appears to be an intrinsic feature of the flaring dynamics, which requires a further 
elaboration of the present solar flaring models. 

It is worth mentioning that we see somewhat intermediate character of the 



reported results with respect to the works by Wheatland and Litvinenko (2002 ), 



Boffetta et al. (1999) and by Lepreti, Carbone, Veltri (2001). The "memory 



revealed in the data rejects models with Poisson statistics. From another per- 
spective, a simple power-law fit very likely fails for the records relevant to a 
specific solar cycle phase. 

The piecewise-constant Poisson model had inspired us to consider solar cycle 
phases separately. We support the argument of this model that solar cycle phases 



9 Recall that memoryless stochastic processes are those with h = const, i.e. the probability of 
an event is independent of neither the waiting time elapsed before it occurs nor the preceding 
event. 
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having quite inhomogeneous statistical properties that are to be considered 
jointly. In other words, the separate consideration of the solar cycles should 
be a cornerstone of a realistic modelling of the solar flaring activity. 

On the other hand, our findings exclude the presence of the timescale of true 
rate constancy (with GOES catalog precision, at least). Even if it exists, it would 
be very small, and may be hard to detect. 

The application of Renewal Theory emphasises the importance of the short 
(energetically small) events in the pdf, whose dynamics had not been pointed 
out in the literature: partially due to removal from the dataset of the B-class 
flares, and partially due to systematic fitting solely the tails of the pdf 's. This 
brings us to problems similar to those that had arisen in the context of coronal 
heating by nanoflares, where, perhaps, the most significant effects are at the 
limit of the instrumental noise. 
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Figure 1. Estimated waiting times versus monthly sunspot number for Ho flares (top) and 
X-ray flares (middle). The sunspot data are provided by the Solar Influences Data Analysis 
Center. The solar maxima (minima) phases that have been used in the analysis are designated 
by thick grey lines (bottom). 
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Figure 2. The pdf of the flare waiting times. The waiting times are estimated from the whole 
GOES Soft X-ray dataset by means of the different wt definitions, which are shown in separate 
panels. The solid line is the pure exponential reference curve with the rate estimated from the 
data being shown by the scatter plot. 



SOLA: wt.tex; 14 September 2012; 0:12; p. 14 



Solar Minimum Phase 




x (min) 

Figure 3. The pdf of the flare waiting times. The waiting times are estimated from the 
GOES Soft X-ray datasct reduced to the years of the solar minimum by means of the different 
wt definitions, which are shown in separate panels. The solid line is the pure exponential 
reference curve with the rate estimated from the data being shown by the scatter plot. 
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Figure 4. The pdf of the flare waiting times. The waiting times are estimated from the GOES 
Soft X-ray dataset reduced to the years of the solar maximum by means of the different 
wt definitions, which are shown in separate panels. The solid line is the pure exponential 
reference curve with the rate estimated from the data being shown by the scatter plot. 
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Figure 5. The estimated flare rate h. The whole set of the GOES Soft X-ray events is used. 
The variation due to wt definition is shown by the scatter plots, which differ from each other in 
the region of short wt . The solid lines are the smoothing applied to emphasise the qualitative 
behaviour. The _B-class flares arc included. 



SOLA: wt.tex; 14 September 2012; 0:12; p. 17 




x (min) 

Figure 6. The scatter plot of the smoothed flare rate h. The whole set of the GOES Soft 
X-ray events is used. The cases without B-class flares are shown for comparison. 
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Figure 7. The smoothed flare rate h for the different solar cycle phases are shown separately. 
The cases with B-class flares excluded are shown by grey lines. See text for discussion. 
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